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Abstract 



We point out that the abelian projection theory of quark confinement is in con- 
flict with certain large-N predictions. According to both large-N and lattice strong- 
coupling arguments, the perimeter law behavior of adjoint Wilson loops at large scales 
is due to charge-screening, and is suppressed relative to the area term by a factor of 
1/N 2 . In the abelian projection theory, however, the perimeter law is due to the 
fact that N — 1 out of iV 2 — 1 adjoint quark degrees of freedom are (abelian) neutral 
and unconfined; the suppression factor relative to the area law is thus only 1/N. We 
study numerically the behavior of Wilson loops and Polyakov lines with insertions 
of (abelian) charge projection operators, in maximal abelian gauge. It appears from 
our data that the forces between abelian charged, and abelian neutral adjoint quarks 
are not significantly different. We also show via the lattice strong-coupling expan- 
sion that, at least at strong couplings, QCD flux tubes attract one another, whereas 
vortices in type II superconductors repel. 
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1 Introduction 



Large-N arguments are a useful way to critique some of the quark confinement mecha- 
nisms that have been proposed over the years. For example, the fluxon mechanism 
|Ij] does not give a string tension for adjoint quarks at any length scale. This was 
shown to be in contradiction with large-N factorization, which predicts that: (i) the 
adjoint string tension a a (from the confinement scale to the charge screening scale L s ) 
is roughly twice the fundamental string tension o>; and (ii) the distance L s where the 
adjoint flux tube breaks, due to charge screening, goes to oo as N — > oo 0. Existing 
Monte Carlo calculations appear to be generally consistent with these large-N predic- 
tions 0. A quite different mechanism is the "dual superconductor" idea, particularly 
as formulated by 't Hooft in abelian-projection gauges Q; this formulation has been 
widely discussed in recent years [5-11]. It is obviously of interest to see whether the 
abelian projection idea is in agreement with large-N predictions. 

At first sight, the abelian projection theory meets the large-N criteria quite well: 
adjoint quarks will indeed have a string tension a a ~ 2o> from the confinement 
scale to some intermediate distance, beyond which the adjoint quarks are unconfined 
|J. However, as discussed in section 3, there is still a profound contradiction between 
large-N and abelian-projection predictions, namely, for an adjoint Wilson loop of area 
A and perimeter V we expect at large-N the leading behavior 



W(C) = N 2 
while for abelian projection we find 

W{C) = N 2 



N 2 



AT J O 
This may not seem like a very serious difference - just a 1/N vs. 1/A 2 suppression 
of the perimeter law term - but in fact it reveals something very fundamental about 
the abelian projection mechanism: In abelian projection gauge, certain adjoint quark 
colors are unconfined not because of charge screening, but rather because they are 
neutral with respect to the residual ^(l)^ -1 gauge invariance, and this in turn leads to 
results that differ with large-N. The disagreement suggests certain numerical tests of 
the abelian projection idea, based on the behavior of Wilson loops and Polyakov lines 
with insertions of various abelian charge projection operators. It is found, in section 
4, that these projection operators do not make much difference so far as string tension 
and screening distance are concerned, and that the force between quarks which are 
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neutral with respect to the abelian projection subgroup is about the same as the force 
between quarks which are charged in the abelian subgroup. These numerical results 
do not favor the abelian projection theory, and add to the negative body of evidence 
already presented in ref. |7|, §, We will comment on the apparently positive results 
presented in ref. ||T0| , |TT . 



We will also comment, in the context of the lattice strong-coupling expansion, on 
the force between QCD flux tubes. Flux tubes in QCD have often been likened to 
Abrikosov vortices in type II superconductors. The force between vortices in type II 
superconductors is repulsive, due to the negative surface energy of the vortices, so 
it is natural to ask whether the force between QCD flux tubes is also repulsive. In 
the next section it will be shown that, at least at strong-couplings, the force between 
QCD flux tubes is actually attractive. 



2 Charge Screening and Flux Coalescence 

Strong-coupling calculations using the Kogut-Susskind Hamiltonian, or based on the 
Heat-Kernel action, as well as Monte-Carlo calculations in the scaling region with 
the Wilson action ||, all give the result that the force between a quark and an 
antiquark, with color charges in representation R of the gauge group SU(N), is 
initially proportional to the quadratic Casimir Cr of the group representation. For 
any representation of SU(N) there is a generalization of the concept of triality in 
SU(3), known as the "N-ality" or, in the mathematical literature, the "class" of the 
representation. Another result of strong-coupling and Monte-Carlo || calculations 
is that at some distance beyond the confinement scale, at a screening distance denoted 
L s , the force between quarks drops abruptly. Beyond the screening distance, the force 
becomes proportional to the smallest Casimir Cl among all representations with the 
same N-ality as representation R. This is known as "charge-screening", and it has 
a simple physical explanation: As the quark-antiquark pair separates and and the 
energy stored in the flux-tube increases, it becomes energetically favorable to pair- 
create gluons. These gluons bind to the quark and the antiquark and, although gluons 
cannot change the N-ality of the quark charge, they can lower the effective charge 
of the quark-gluon state to a representation L of the same N-ality, with the lowest 
possible Casimir Cl- A particular example is the case of having a pair of quarks in 
the adjoint representation. The adjoint representation has N-ality = 0, the same as 
a singlet. A pair of gluons, binding to each quark, can reduce the charge to a singlet, 
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and therefore just break the flux-tube between the adjoint quarks. The mechanism 
is exactly the same as in QCD with dynamical quarks, where a flux-tube is broken 
by quark-antiquark production. For this reason, in pure QCD, we expect an adjoint 
Wilson loop to have an area law from the confinement scale up to a distance L s , 
after which it is screened to perimeter law. This is indeed what is seen in numerical 
experiments [[|. 

An important point, especially for the purposes of this article, is that diagra- 
matically the charge- screening process is always suppressed by a factor of 1/N 2 in 
pure QCD; this is true whether the process is represented by lattice strong-coupling 
diagrams, as shown in Fig. 1, or by high-order Feynman diagrams. Feynman dia- 
grams can be classified according to their associated powers of 1/N 2 , and diagrams 
of leading order - the planar diagrams - satisfy the factorization property 

<(TrA)(TrB)>=<TrA><TrB> (3) 

For this reason a a = 2<tf in the large-N limit; it is also the reason why charge 
screening must be a non-planar process. Non-planar diagrams, whether Feynman or 
lattice strong-coupling, are suppressed (in pure QCD) by a factor of at least 1/N 2 
relative to the planar diagrams. Therefore, we conclude on these very general grounds 
that for an adjoint Wilson loop 



W A (C) = N 2 



(4) 



In strong-coupling one can easily calculate where the perimeter behavior of the second 
term takes over from the area behavior; for square L x L loops this occurs for 



L > 



1/2 

(5) 



_ CTp 

which diverges as iV — > oo. 

An effect which is very closely related to charge screening is the phenomenon of 
flux coalescence. This effect is relevant to the question of whether the force between 
QCD flux tubes is repulsive, as it is for Abrikosov vortices in type II superconductors, 
or attractive. 

Let us consider two quark-antiquark pairs with the axes through each pair parallel; 
these are represented by the two parallel Rx T Wilson loops shown in Fig. 2, which 
are separated by a distance L » R, and which have the same orientation. Denote the 
expectation value of the product of these loops by D(R,T, L). Then the interaction 
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energy between the loops, which represents the interaction of the flux tubes between 
the quark- ant iquark pairs, is just 



d 

V(L) = - lim — In 



D(R,T,L) 
D(R, T,oo) 



and 



D(R, T, L) = D(R, T, oo) + RT{L + l)(A a - A b ) 
where, to leading order 

D(R,T, oo) = N 2 e- 2apRT 



(6) 



(7) 



(8) 



The leading strong-coupling contributions to A a and Aj, are indicated schematically 
in Fig. 3. 

We will calculate V(L) using the heat-kernel lattice action 



v 



(9) 



where the product extends over all oriented plaquettes p and the summation runs over 
all inequivalent irreducible representations r of dimension d r . Xr(U p ) is the character 
of the group element U p G SU (N) in the representation r, and C r is the eigenvalue of 
the quadratic Casimir operator in this representation. For quarks in the fundamental 
(defining) representation, the plaquette with character xt(U P ) in the middle of the 
tube in Fig. 2a can only be in the symmetric (r = s) or antisymmetric (r = a) 
representations formed from the product of two fundamental defining representations. 
The interaction potential is then found to be 



V(L) = —(L + l)e~ 4r7FL R 



N-l A + N. N+l 1 - N ' 



(10) 



2N ~""" L N 2 (3 ' 2N N 2 (3 
Expanding the exponents in a Taylor series, V(L) to leading order in 1/N is given by 

1 1 



V{L) 



2f3 2 N 2 



R(L + l)e~ 4aFL 
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which is an attractive potential at large L » Flux tubes therefore tend to 
attract one another, by a process which can be viewed as glueball exchange. 

For separations L « R flux tubes do not only attract, they will actually tend to 
coalesce into a single flux tube. In this case 

D(R, T, L) = D(R, T, oo) + D C (R, T, L) (12) 

where the leading contribution to D c is shown in Fig. 4. In this figure there is a single 
sheet of plaquettes, of area R x T, in a fundamental representation of N-ality N-2. 
For R « NL the disconnected diagrams dominate, and the energy of the system is 
simply the sum of the two disconnected flux tubes 

N 2 - 1 

E = 2a F R = -j^p-R (13) 

whereas, if R » NL, it is energetically favorable for the two flux tubes to coalesce, 
as represented in Fig. 4, and the energy is 

(N-2)(N + 1) N 2 - 1 
E c = a a R + 2a F L = [ - j$^R + -j^-L (14) 

It is clear that E c < Eq, i.e. coalescence is favored, for R > NL. 

The calculation above can be easily generalized to the case of n quark- ant iquark 
pairs, with axes nearby and aligned parallel to one another. Once again, for quark 
separations L much less than quark-antiquark separations R, it is energetically fa- 
vorable for the n flux tubes between each quark-antiquark pair to coalesce into a 
single tube, as shown in Fig. 5. Its string tension is determined by the eigenvalue of 
the quadratic Casimir operator of the lowest dimensional representation r of N-ality 
(N — n)modN, which is again a fundamental representation of dimension d r = (^j . 
It is clear that coalescence may result in considerable lowering of the energy of the 
system, compared to the configuration of n separate flux tubes between each qq pair; 
the ratio between the energies of these two different configurations being, for R » L, 

ol C L , n-1 



nap tiCf N — 1 

For n=2 in SU(3), this ratio is just r = |. 

The phenomenon of flux-tube coalesence implies that the force between nearby 

QCD flux tubes is attractive, at least at strong couplings, in contrast to the situation 

3 In fact, this potential is also attractive at small L, given of — Cf/{N(3) > 1, which is the case 
at strong-couplings. But at L « R, flux tube coalescence, discussed below, is the dominant effect. 
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in type-II superconductivity. If strong-coupling calculations are any guide, and if 
confinement in QCD is analogous to (dual) superconductivity, then the analogy would 
presumably be to type I superconductivity. More importantly, both color charge 
screening and the attractive force between flux tubes are due to non-planar, 1/N 2 
suppressed processes. We will now consider whether the abelian projection theory is 
consistent with this suppression factor. 

3 The Abelian Projection Theory 

The abelian projection theory of quark confinement was put forward in ref. 0] . In this 
theory a gauge-fixing condition is first chosen to break the SU(N) symmetry down 
to the Cartan subgroup [/(l)^ -1 . Monopoles are then identified with singularities 
in the gauge-fixing condition, and condensation of these monopoles is invoked to 
explain confinement of particles charged with respect to the ^/(l)^ 1 subgroup. The 
confinement mechanism is analogous that of compact QED, with gluons associated 
with the Cartan subalgebra playing the role of the photon field which forms the flux 
tube. These "diagonal" gluons are uncharged with respect to the residual [/(l)^ 1 
symmetry; all other gluons are charged with respect to the residual symmetry, and, 
like quarks, are confined by flux tubes. 

Consider the force between two quarks with charges in adjoint representation, in 
the abelian projection theory. For simplicity (and because it makes little difference 
in the large-N limit), take the group to be U(N) rather than SU(N), so the abelian 
subgroup is U(1) N . Let gi denote an abelian charge of magnitude g in the i-th U(l) 
subgroup. A quark in the fundamental (defining) representation has iV color degrees 
of freedom, and from the abelian-projection point of view, each of these degrees of 
freedom corresponds to a quark Q % F with a different abelian charge g^, 1 < i < N. 
Likewise, the N 2 degrees of freedom of a quark in the adjoint representation can be 
be grouped according to abelian charge: there are N(N — 1) quarks Q 1 ^, i ^ j with 
charge (g^ —Qj), and there are iV quarks Q\ which are neutral with respect to the 
U(1) N subgroup. 

Since the flux-tube between static adjoint quarks is neutral with respect to the 
abelian subgroup in the abelian projection picture (flux tubes are formed from the 
"photon" fields), the charges of the quarks themselves must be correllated in order 
to have a neutral composite state. The quark-content of the composite state, in an 
abelian-projection gauge, must then be 
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TrQ[ Xl ]Q[x 2 ] =J^QV[xi]<y i [x 2 ]+^Q a [x 1 ]Q«[x 2 ) (16) 

ijtj i 

Consider the N(N — 1) contributions to the first sum. Each contribution represents 
two quarks with equal and opposite charges (gi, —gj) and (— gi, gj) respectively. Since 
each quark is charged in two different abelian subgroups, there will be two flux tubes 
between them, each with string tension op. Thus the net string tension for these 
quarks is a a = 2a>, in agreement with large-N. The N quarks of the second sum, 
however, are neutral with respect to the confining abelian U(1) N subgroup; these 
neutral quarks are unconfined. For an adjoint Wilson loop of area A and perimeter 
V, we therefore expect 

W(C) = N(N - 1) exp[- 2a F A - fiP] + Nex V [-^iV] (17) 

for the following reason: Given that the abelian "charged" gluons are confined and 
the flux tube is neutral, the quarks cannot exchange charge between them, and a 
quark can only change its abelian charge by emission and reabsorbtion of a charged 
gluon. Even allowing for such virtual processes, the abelian charge of each quark + 
gluon-cloud cannot change. Thus the color sum in a Wilson loop can be expressed 
as a sum of loops, each representing a particular (gi, —gj) charge running around the 
loop. The N(N — 1) charged adjoint quarks (i ^ j) contribute to the first term in 
(|TT|), and the N neutral, unconfined quarks contribute to the second term. 

The expression ([17]) above for the adjoint Wilson loop is in disagreement with the 
large-N prediction. In the large-N analysis, all quark colors are on the same footing, 
and the perimeter term comes about through charge- screening, and not because any 
subset of quark charges is oblivious to the confining force. More importantly, accord- 
ing to large-N, the coefficient of the perimeter term is 0(1), rather than 0(N) as in 
the abelian projection theory. 

To make the same point in a slightly different way, consider the trace of an adjoint 
Wilson loop in SU(N) 



Tr A 



Pexp[i f A a L a ] = J2 Tr A \m){m\P exp[i f A a L a ] 
J J m I j 

= ^2(m\P exp[i j> A a L a ]\m) 



where L a are the group generators and sum is over all members \m) of the multiplet. 
The confining force can only come from coupling of the adjoint quarks to the gluons 
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in the Cartan subalgebra, with generators denoted Hi (i = 1, N). Therefore, in an 
abelian projection gauge such as the maximal abelian gauge, following the reasoning 
ofref. 0, 



< Tr f 



Pexp[i j> A a L a 



> 



<Tr A 

N 2 -l 



exp[i I A 1 Hi 
E < ( m \ exp[i <j) A l Hj\\m) > 

n=l 

E <exp[^^ m) ]> 



(19) 



where the are eigenvalues 



m=l 



Hi\m) = |ra) 



(20) 



Then, making use of the fact that the multiplicity of zero-weight (all \ = 0) states 
in the adjoint representation is N — 1, we have 



+ (N-l) 



(21) 



N(N-l) 

W{C) « £ 

m = 1 
(non-zero 
weight states) 

Ignoring the coupling of adjoint quarks to gluons which are not in the Cartan subal- 
gebra, which was the approximation made in eq. (0), has the effect of dropping self- 
energy (perimeter law) contributions to the zero-weight and underestimating them 
for the non-zero weight states. If these contributions are included perturbatively, the 
final answer would be 

N(N-l) 

W(C)n £ e -^A-nV + (jv-l) e -^ (22) 

m = 1 
(non-zero 
weight states) 

Again, the deconfined term is O(N), rather than 0(1) as expected by large-N argu- 
ments. 
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4 Numerical Tests 



The fact that the abelian projection theory disagrees with large-N predictions should 
not be too surprising. As noted above, the deconfinement of adjoint loops according 
to the abelian projection is simply due to the fact that N — 1 out of N 2 — 1 adjoint 
quarks are neutral with respect to the [/(l)^ 1 subgroup, and are therefore unconfmed 
in the abelian projection theory. This gives a suppression factor of only 1/N to the 
unconfmed contribution. Charge-screening, which involves pair-creating gluons that 
bind to the adjoint quarks, is a completely different mechanism, and just gives the 
usual non-planar suppression factor of 1/N 2 . This disagreement doesn't necessarily 
mean that the abelian projection theory is wrong; it could be the large-N arguments 
that have somehow gone astray. But such a clear difference in the two approaches 
does suggest a simple numerical test. 

Let us consider the case of N = 2, with H = L3 to be the generator of the 
Cartan-subalgebra. Then eq. ( |19D becomes 

W{C) « ]T < exp[tm j> A 3 } > 

m=— 1 

= 2e aWA + l (23) 
or, correcting for the self-energy (perimeter) contributions 

W(C) = 2e aWA -» v + e-* v (24) 

In other words, the abelian-projection prediction is that, in an abelian projection 
gauge 



< (l|Pexp[i j A]\l) > = < (-l\Pexp[ij> A}\ - 1) >= e -^ w ^-^ 

< (0|Pexp[z^A]|0) > = e~^ v (25) 

The inner product < (0|...|0) > can be thought of as creation of a quark- ant iquark 
pair with zero abelian charge, which run around the loop. The only way to change the 
abelian charge of the pair would be for the quarks to exchange a gluon with non-zero 
abelian charge. However, this process should be suppressed at the confinement scale, 
since the charged gluons are presumably confined, and therefore only contribute to 
the self-energy of each quark. The string tension of the < (0|..|0) > term should 



11 



therefore be zero in an abelian-projection gauge; there should be no area suppression 
whatever, in any range of quark separations. This is an easy prediction to test. 

Let U (C) represent a product of link variables along the path C. In SU (2) this 
can always be expressed as 

U(C) = a l + ia k a k (26) 
where the a k are the Pauli matrices. Then, in the adjoint representation 



U ++ {C) = {l\U A {C)\l) = a 2 -al + 2ta a 3 
C/__(C) = (-1\U A (C)\ -1) = a 2 -al-2ia a 3 

Uoo(C) = (0\U A (C)\0) = 2(al + a 2 3 )-l (27) 



Define the charged and neutral Wilson loops 



W C (C) = < U ++ (C) >=< C/__(C) > 

W (C) = <U 00 (C)> (28) 

together with the corresponding Creutz ratios Xc[R,T] and Xo[R, T\. As discussed 
above, the abelian projection prediction is that Xo[R,T] = (or, at least, that 
Xo << Xc) in an abelian projection gauge. 

We have computed these Creutz ratios by lattice Monte Carlo, in D=3 dimensions 
with a Wilson action at (5 — 5, which is just inside the D=3 scaling region. The 
charged and neutral loops were evaluated in maximal abelian gauge. The results 
obtained after 153000 update iterations, for x [i?, -R], Xc[R,R], and for x[Ri R] (the 
Creutz ratio of gauge-invariant adjoint j '• = 1 loops), are shown in Fig. 6. Creutz 
ratios for the fundamental (j = |) loop xf(R,R) are also displayed in Fig. 6. It 
can be seen that Creutz ratios in the j — \ and j = 1 representations are in the 
proportion predicted by large-N, which is a ratio of Casimirs 

x[R,R\ ^ (29) 



Xf[R,R] 3 

It can also be seen that Xo is n °t zero, and shows no tendency to go zero faster 
than Xc- In fact, although xo is smaller than Xc, the difference is only about 10%; this 
difference does not appear to grow with loop size. Thus, while there may perhaps be 
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some effect attributable to abelian monopoles, the effect seems rather small; certainly 
it is not sufficient to explain confinement in this region. This result seems entirely 
consistent with the results in ref. ||. 

We have also studied the correlation of Polyakov lines in D=4 dimensions and 
maximal abelian gauge 

P (R) = <U Q0 [L 1 )U 00 [L 2 ]> 

P C (R) = <U ++ [L 1 ]U—[L 2 ]> (30) 

where L 12 are parallel Polyakov loops, 2 lattice spacings in length, separated by R 
lattice spacings in the spatial hyperplane. Because of the small time extension, we 
must work at a strong coupling - in our case (3 = 1.8 - to avoid the deconfinement 
transition. Although this is at strong-coupling, it is at least close to the strong- 
to-weak coupling transition point, and we can see if there is any tendency for the 
"neutral" correllations Po(R) to behave differently from the "charged" correlations 
P C {R)- The data is shown in Fig. 7. Charge screening sets in at about 2 lattice 
spacings, and there does not appear to be any difference in the behavior of Po and 

Pc 

For the numerical simulations, we used an heat-bath algorithm, running on a 
16 3 x 2 lattice at (3 = 1.80. The results were obtained averaging over 1000 config- 
urations, after 20000 sweeps of thermalization. The continous lines in Fig. 7 are the 
results of fit for the thermodynamical mixing of two interaction channels as described 
in ref. ]Tj|. For the first Boltzmann factor we used the lattice version of a Coulomb 
plus linear potential plus self energy and via the second term we took the screening 
of charges into account. 

The conclusion of this numerical work is that there seems to be no appreciable 
difference in the forces between abelian charged and abelian neutral adjoint quarks, 
at the couplings and separations we have investigated. This evidence does not favor 
the abelian projection theory. 

It may be appropriate, at this point, to comment on certain other Monte Carlo 
investigations of the abelian projection theory. Numerical work on this problem was 
initiated by Kronfeld et. al |5|, who found a drop in the monopole density at the 
deconfinement phase transition, in "maximal abelian gauge" defined as the gauge 
which maximizes the quantity 

Q = EE^ht/,w#;wi (si) 

X [1=1 
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More recent investigations by Del Debbio et. al. have shown, however, that the 
definition of monopole density is plagued by lattice artifacts, and is not at all a good 
order parameter for confinement |7|]. For example, it is found the monopole density 
neither shows correllation with the string tension with cooling, nor a drop across 
the deconfinement phase transition. These problems may be alleviated by a more 
appropriate definition of the monopole creation operator, but we will not pursue that 
issue here. 

An alternate line of investigation is to see if the electric flux inside the flux tube 
is dominated by the Cartan subalgebra; e.g. in SU(2), with the choice of maximal 
abelian gauge above, one checks to see if the field-strength is proportional to 03. 
Let Upu be the plaquette variable. One defines the field strength on the lattice as 
F" = Tr\U^ v — U^a" 1 ]/ '4i, as well as the quantities 

< D a TrW > 

J a = J TI - <D a > 32 
< TrW > v ; 

where W is a Wilson loop in the pv plane; and D a = J2{F^(x)) 2 /np, where the sum 

is over central plaquettes in the minimal area bounded by the loop, furthest from the 

boundary, np is the number of central plaquettes. Let J = J 1 + J 2 + J 3 , and define 

the ratio 

P = y (33) 

If p ~ 1 this would tend to confirm the abelian projection theory, while p ~ 1/3 
would tend to refute it. The result found in |8|, |9| was p ~ 1/3. 



On the other hand, Hioki et. al. |T0[ have argued that this ratio should be defined 
using abelian loops in the numerator. An abelian link is defined as the diagonal part 
of the link variable, rescaled to restore unitarity. An abelian loop is a loop constructed 
from abelian links. Let be a 1 x 1 abelian loop, and define 



Ja 

PA = T 



< — E fLTrW > 
< TrW > 

U» = (u, v -u; v )/2i (34) 



j a = : ^„ TT7 : < jpv > 



With these definitions, Hioki et. al. find pa actually greater than 1 (up to pa ~ 1.6, 
in fact). 
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We regard conclusions based on the enhancement of abelian loops, both in ref. 
10| and also in [|llj as very misleading for the following reason: Maximal abelian 



gauge simply makes links as diagonal as possible; the enhancement of abelian loops 



(termed " abelian dominance" ) which was found in IJlTJ is a simple consequence of this 
fact. It is therefore crucial to choose observables whose behavior will really test the 
abelian projection theory, rather than just display this particular aspect of the gauge 
condition. The behavior of abelian Wilson loops, and the pa quantity defined above, 
do not meet that criterion. 

The following calculation will illustrate the point. Instead of gauge-fixing to max- 
imal abelian gauge, let us fix to another, "x-y" maximal abelian gauge, introduced in 
ref. ||, which is defined to maximize the quantity 

Q = Y.il Tr \°zU,(x)azUl(x)] (35) 

x ^=1 

This gauge forces links in the x and y directions only to be as diagonal as possible. 
Since there is no requirement, in the abelian projection theory, that the gauge-fixing 
condition must be spherically symmetric, this gauge should be as good as the usual 
maximal abelian gauge. We now compute pa by Monte Carlo separately for loops 
oriented in the x-y, and z-t planes. The results have only been computed for rather 
small loops, but we feel they already show how things go. The loops are from R=l to 
R=4 lattice spacings wide, and T=4 lattice spacings high. Half of the plaquettes in the 
minimal area of the loop were used for calculating J a] these are the 2 x R plaquettes 
which are one lattice spacing away from loop boundaries in the T direction. 

The results of this calculation, performed in D=4 dimensions at /3 — 2.4, are 
shown in Fig. 8. If one accepts the proposition that pa ~ 1 is evidence for the abelian 
projection theory, then it would seem from this data that monopoles are responsible 
for confinement in the x-y plane, but not in the z-t plane. That conclusion, we feel, 
is nonsense. A much more reasonable explanation is that, because links are nearly 
diagonal in the x-y plane, loops which are built from the diagonal part of the links 
pretty well approximate the full Wilson loop. Conversely, in the z-t plane, the results 
look much as if there were no gauge-fixing at all. 

As an additional check, we have also computed pa for a small R = 2, T = 4 loop 
at f3 — 2.8, in ordinary maximal abelian gauge. At this value of f5, for such a small 
loop, there should be no flux tube formed; yet we find pa = 0.88. This is another 
indication that the large value of pa is a gauge effect, which has nothing to do with 
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flux tube formation. 

This calculation illustrates the fact that, in maximal abelian gauge, one must 
ensure that the observables chosen are relatively insensitive to the diagonality of 
links which is enforced by the gauge. This is the case for the J a observables of eq. 
(j32|), whose sum is directly related to the energy density in the flux tube; it is not 
the case for the J a observable of eq. (|51D, which is not related directly to the energy 
density. 



5 Conclusions 

A prediction of the abelian projection theory of confinement is that, in an abelian 
projection gauge, adjoint quark colors which are neutral with respect to the remnant 
C/(l) Ar_1 symmetry are oblivious to the confining force. This prediction turns out to 
conflict with large-N arguments, and can be tested by looking for the absence of an 
area-law term, over any length scale, in Wilson loops and Polyakov lines with ap- 
propriate insertions of abelian-neutral projection operators. We have found, instead, 
that insertion of abelian neutral (and abelian charged) projection operators has very 
little effect on the value of the string tension extracted from the loop. Assuming that 
the QCD flux tube in abelian projection gauge is abelian neutral (since it is supposed 
to be formed by the " photon" fields) , this means that there is no significant difference 
in the forces between abelian neutral adjoint quarks, and between abelian charged 
adjoint quarks. Thus our data is consistent with large-N expectations, and in conflict 
with the abelian-projection theory, which holds that the confining force is sensitive 
mainly to the U{l) N ~ l charge. 

We have also examined the claim of "abelian dominance" in maximal abelian 



gauge, found in [pi], |TT|. We believe that our data in an "x-y maximal abelian gauge" 
(together with previous work along these lines in ref. [§]), demonstrates that this 
enhancement of abelian loops is purely a gauge effect, with no relevance at all to the 
physics of confinement. 

Finally, we have calculated the force between QCD flux tubes in the lattice strong- 
coupling expansion, and find that this force is attractive. If the strong-coupling result 
survives in the continuum theory (and there is reason, based on the general notion of 
flux coalescence, to think that it might), then QCD flux tubes are not analogous to 
Abrikosov vortices in type II superconductors, which tend to repel one another. 

In general, theories of confinement which rely on analogies to abelian gauge theo- 
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ries will tend to identify a small subset of the degrees of freedom, e.g. associated with 
the Zjsr center or the Cartan subgroup of the full gauge group, as being especially 
important for quark confinement. Such a subset becomes negligibly small compared 
to the total number of degrees of freedom in the iV — > oo limit, and it is therefore 
not surprising that such theories will somewhere contradict results based on large-N 
counting arguments. One theory of confinement which is based rather explicitly on 
the large-N picture is the "gluon-chain" model of flux-tube formation, advocated by 



one of us in ref. [fL4}l . This model is consistent with all large-N predictions, and also 
has some numerical support |I5|. It is relevant here as an example of a confinement 
mechanism in non-abelian gauge theories which has no abelian counterpart; we believe 
that this must be true of any mechanism which is consistent with large-N results. 
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Figure Captions 



Fig. 1 Leading strong-coupling contributions to the adoint Wilson loop, (a) the 
0(N 2 ) area contribution; (b) the 0(1) perimeter contribution. 

Fig. 2 Two parallel quark- ant iquark flux tubes realized by two parallel Wilson loops 
of the same orientation. 

Fig. 3 Diagrams used in computing flux-tube interactions (eq. (0)): (a) term A a , (a 
plaquette in the middle of a tube); (b) term Af,. It is necessary to sum over 
the position of the tube in the sheet, and over the position of the "middle" 
plaquette along the tube. 

Fig. 4 Leading contribution for parallel Wilson loops at L << R. The sheet in the 
middle (bounded by the heavy solid line) is in the fundamental representation 
of N-ality N-2. 

Fig. 5 Coalescence of n quark-antiquark flux tubes into a single flux tube of N-ality 
(N - n)modiV. 

Fig. 6 Creutz ratios of adjoint (J = l)Wilson loops. Stars are the usual % ratios, 
diamonds are ratios xo of loops with abelian-neutral (m = 0) projection op- 
erators; vertical crosses are ratios Xc of loops with abelian-charged (m = ±1) 
projection operators, evaluated in maximal abelian gauge. The errors on loops 
at R = 1, 2, 3a are negligible, errors at R = 4a are about ±10%. Also shown 
(by diagonal crosses) is xf, the usual Creutz ratio in the fundamental (J = |) 
representation. Simulation is in D=3 dimensions at (3 — 5. 

Fig. 7 The potential V(R) extracted from Polyakov lines of extension N t = 2 lat- 
tice spacings. The potentials for "neutral" (diamonds) and "charged" (vertical 
crosses) sources in maximal abelian gauge of SU(2) are compared with the un- 
gauged potentials (stars) in the adjoint representation. For comparison we show 
also the potential between doublet sources. The statistical errors start to get 
larger than the symbols in the screening region and can be estimated from the 
fluctuations with the distance R. The full lines are the results of fits with a sum 
of two Boltzmann factors. 
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8 pa in the "x-y maximal abelian" gauge. Crosses represent pa extracted from 
loops in the x-y plane; diamonds represent pa extracted from loops in the z-t 
plane. Squares are values of pa taken with no gauge fixing. 
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